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I .  Introduct ion 


The  inductive  acceleration  of  electrons  by  adiabatic  magnetic  com¬ 
pression  of  electron  rings  has  been  demonstrated  in  several  lanorator  ies . 1 
For  weak  rings,  the  observed  increase  of  the  mean  electron  energy  and  the 
shrinkage  of  the  major  and  minor  radii  of  the  ring  are  consistent  with 
the  constancy  of  the  magnetic  moment  j±[=  P^/2mQB,  where  P  is  the  single 
particle  momentum  perpendicular  to  the  magnetic  field  and  B  the  magnetic 
field].  A  two-hundred  fold  increase  of  the  electron  energy  has  been  ob¬ 
served  by  Kapetanakos2  et  al.,  when  the  time  varying  magnetic  field  in¬ 
creased  from  its  initial  value  of  70  G  to  15  kG. 

The  inductive  acceleration  of  ions  by  magnetic  compression  of  ion 
rings  is,  in  principle,  very  similar  to  that  of  the  electrons.  However, 
in  practice,  the  ion  ring  compression  is  complicated  by  the  presence  of 
space  charge  neutralizing  electrons.  These  electrons  can  generate  an 
azimuthal  current  either  by  electron-ion  drag5*6,  radial  electric  fields 
or  pressure  gradients.  Thus,  depending  on  the  conditions  of  the  system, 
at  least  partial  cancellation  of  the  current  that  is  carried  by  the  ions 
is  likely,  in  particular  near  the  magnetic  field  null.  In  some  application 
partial  cancellation  of  the  ion  current  during  magnetic  compression  is 
highly  desirable  because  most  of  the  externally  supplied  energy  goes  into 
useful  kinetic  energy  and  not  into  thermal  and  self-magnetic  field  energy 
of  the  ring. 

During  the  last  few  years,  several  models7-10  have  been  developed  in 

relation  to  the  magnetic  compression  of  ion  rings.  The  predictions  of 

these  theoretical  models  are  sensitive  to  their  basic  assumptions.  In  all 

these  models  the  dynamics  of  the  neutralizing  electrons  either  have  been 
Note  Manuscript  submitted  August  1,  1979. 
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neglected  or  treated  incompletely. 

In  a  recent  publication  Sudan11  has  suggested  the  use  of  magnetically 
compressed  rings  for  pellet  fusion.  In  Sudan's  scheme  the  image  currents 
on  the  walls  of  the  tube  that  surrounds  the  ring  provide  the  radial  equi¬ 
librium  during  the  propagation  of  the  ring  from  the  compression  region  to 
the  target.  The  guide  tube  is  destroyed  and  must  be  replaced  in  each  shot. 

In  this  report,  we  propose  an  alternative  scheme  for  the  generation 
of  high  energy,  high  current  ion  pulses  that  also  is  based  on  the  magnetic 
compression  of  ion  rings. 

The  basic  difference  between  the  present  scheme  and  that  of  Sudan 
is  on  the  extraction  and  the  propagation  of  the  ion  pulse  after  compression. 
The  proposed  scheme  does  not  require  either  an  external  magnetic  field  or 
a  tube  for  guiding  the  ion  pulse  from  the  compressor  to  the  target. 
Therefore,  it  may  be  in  particular  useful  in  those  applications  that  require 
an  appreciable  separation  between  the  compressor  or  accelerator  and  the 
target.  In  addition,  the  present  approach  assumes  that  the  compression 
of  the  ion  ring  is  performed  in  a  suitable  environment  that  allows  ap¬ 
preciable  cancellation  of  the  ion  current. 

The  proposed  scheme  for  the  generation  of  intense,  high  energy  ion 
pulses  by  magnetic  compression  of  ion  rings  is  shown  schematically  in 
Fig.  1.  A  hollow,  thin  ion  beam  of  energy  about  2  MeV  generated  by  a 
low  inductance,  coaxial  reflex  tetrode12  is  passed  through  a  magnetic 
cusp.  As  a  result  of  the  q(vz  x  Br)  force,  the  ions  of  the  beam  obtain 
azimuthal  velocity  and  start  to  rotate.  The  rotational  velocity  (vq)  of 
the  ions  is  further  enhanced,  at  the  expense  of  their  translational  velo¬ 
city,  by  a  static,  converging  magnetic  field.  The  maximum  value  of 


2 


the  compressing  field  is  such  that  the  protons  located  at  the  outer  edge 

of  the  beam  arrive  at  B  with  zero  translational  velocity.13  The  ion 

max 

ring  is  formed  by  trapping  the  50-70  nsec  duration  ion  pulse  inside  a 
magnetic  mirror  with  the  aid  of  a  gate  field. 

The  rotational  energy  of  the  r'ng  is  enhanced  by  increasing  in  time 
the  confining  magnetic  field.  For  adiabatic  compression,  an  appreciable 
saving  in  magnetic  energy  is  obtained  by  using  imploding  liners.14  After 
compression,  the  ring  is  extracted  out  of  the  confining  field  by  opening 
the  far  mirror  peak.  Initially,  the  ring  is  allowed  to  expand  adiabatical ly 
in  a  spatially  decreasing  field.  When  the  ratio  v„/v  ,  where  v„,v 

i  i 

are  the  velocities  of  the  ring  parallel  and  perpendicular  to  the  field 
lines  respectively,  acquires  the  desired  value,  the  ring  is  passed  through 
a  sharp  half  cusp  that  further  increases  v„  at  the  expense  of  v  . 

Although  the  radius  of  the  beam  remains  virtually  unchanged  as  it  passes 


through  the  sharp  half  cusp,  the  conservation  of  the  canonical  angular 
r  9rA„-i 


momentum  P  =  mrv„  +  — 
9  c 


requires  a  rapid  expansion  of  the  beam  when 


Ag(r)  =  0  .  However,  for  intense  rotating  beams  Ag(r)  4  0  on  the  right 

of  the  half  cusp  because  Ag(r)  =  AgXt(r)  +  Age'^(r)  and  Age'^(r)/  0  ,  al¬ 
though  the  externally  applied  vector  potential  AgXt(r)  =  0  .  The  existence 
of  equilibrium  for  a  hollow,  rotating  ion  pulse  that  propagates  in  the 
A„  (r)  =  0  region  is  shown  in  Section  IV,  using  the  Vlasov-Maxwel 1 


equat i on . 


The  organization  of  this  report  is  as  follows:  The  magnetic  compression 
of  the  ring  is  discussed  very  briefly  in  Section  II.  The  transmission  of 
the  ring  through  the  cusp  and  the  target  location  for  the  particular  case 
of  pellet  fusion  is  presented  in  Section  III.  The  equilibrium  of  the 
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hollow  beam  in  the  absence  of  an  externally  applied  field  with  and  without 
a  conducting  tube  is  discussed  in  Section  IV.  Finally,  a  brief  summary 
of  the  results  is  given  in  Section  V. 

II.  Magnetic  Compression  of  the  Ring 

Efficient  transfer  of  the  externally  supplied  magnetic  energy  into 
particle  kinetic  energy  during  magnetic  compression  of  an  ion  ring  requires 
that  the  self  magnetic  field  to  be  appreciably  smaller  than  the  applied 
magnetic  field.  Under  these  conditions,  the  magnetic  field  lines  are  open 
and  thus  only  small  electric  fields  can  be  sustained.  If  the  Budker's 
parameter  V  (=  NRQ,  where  N  is  the  number  of  ions  per  unit  length  and  PD  is 
the  ion  classical  radius)  of  the  ring  is  not  appreciably  smaller  than  unity, 
the  self  field  can  be  kept  small  only  if  the  electron  return  current  part¬ 
ially  cancels  the  ring  current.  Such  an  electron  current  can  be  driven 
by  a  radial  electric  field,  a  plasma  pressure  gradient  along  the  minor 
radius  of  the  ring  or  collisions.  Among  the  three,  plasma  pressure  grad¬ 
ients  appear  to  be  the  most  appropriate  in  the  present  case.  The  radial 
force  required  for  the  existence  of  equilibrium  is  mainly  provided  by  the 
magnetic  field. 

If  the  self  magnetic  field  of  a  space  charge  neutral  ring  is  ap¬ 
preciably  smaller  than  the  external  field,  the  kinetic  energy  of  the  ions 
E(t)  and  the  particle  current  I ( t )  can  be  estimated  from 

E  (t)  =  2E(o)[B(t)/B(o)][Y(t)  +  1]-»,  (0 

and 

Kt)  =  lo(o)[B(t)/B(o)]Y-1(t),  (2) 
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where  E(0),  I (o)  and  B(o)  are  the  initial  values  of  energy,  particle 
current  and  magnetic  field  respectively,  B(t)  is  the  value  of  magnetic 
field  at  time  t  and  y ( t )  the  relativistic  factor. 

The  technique  of  generating  high  magnetic  fields  by  magnetic  flux 
compression  is  under  investigation  for  several  years.  At  NRL,  large 
initial  diameter  (—28  cm),  7“cm  long  seamless  liners  of  annealed  aluminum 
have  successfully  imploded  and  peak  magnetic  fields  in  excess  of  1.3  M 
gauss  have  been  obtained.  In  addition,  experiments  have  shown  that 
rotating,  hollow  cylinders  of  liquid  liners  ( 22 %  Na,  78%  K)  can  be  stably 
imploded  on  trapped  magnetic  flux. 

Presently,  experiments  are  underway  at  NRL  to  form  a  strong  proton 
ring.  The  anticipated  parameters  of  such  a  ring  are  shown  in  Table  1. 

The  values  of  the  various  parameters  after  compression,  are  also  shown 
in  the  table.  In  an  experiment,  the  values  of  the  various  parameters  at 
the  peak  of  the  compression  will  not  be  much  different  than  those  given  in 
Table  I  when  the  self  magnetic  field  is  small  in  comparison  with  the  ex¬ 
ternally  applied  field. 

111.  Extraction  of  the  Ring 

After  the  ring  compression  is  over,  the  far  mirror  is  switched  off 
and  the  ring  is  allowed  to  propagate  adiabatical ly  in  a  slowly  decreasing 
magnetic  field.  The  velocity  of  the  ring  parallel  to  the  magnetic  field 
lines  at  point  1  (cusp  entrance)  is  given  by 

v„i  *  v10(l-B1/Bm)1/2,  (3) 

where  vH  is  the  azimuthal  ring  velocity  at  the  end  of  the  compression. 
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B  is  the  middle  plane  value  of  the  compressing  field  and  B,  is  the  magnetic 
m 

field  at  the  entrance  of  the  cusp.  The  ratio  B,/B  is  mainly  dictated  by 
the  maximum  ion  pulse  radius  that  can  be  tolerated  in  a  particular  appli¬ 
cation.  For  example,  in  pellet  irradiation  experiments  with  the  parameters 
of  Table  1,  Bj/B^  =  1/2  appears  to  be  appropriate.  In  order  to  further 
increase  v  without  an  appreciable  radial  expansion,  the  ion  pulse  is  passed 
through  a  half  cusp  that  has  a  transition  width  6j  that  is  smaller  than  the 
average  radius  of  the  gyrating  ions. 

It  is  assumed  that  the  cusp  is  described  by  the  vector  potential 

Aq  =  rzBj/Ufij)  ,  *  6,  <  z  <  0  ,  (M 

in  a  cylindrical  coordinates  system  with  its  origin  as  shown  in  Fig.  1.  For 
a  charged  particle  of  mass  m  and  charge  q  in  an  axisymmetric  field  the  non- 
relativistic  Hamiltonian  is 

H  =  (m/2)  (v£  +  v0  +  v*)  +  qif>(r)  ,  (5) 

where  <)>(r)  is  the  electrostatic  potential  and 

q  rAfl 

ve  =  (pe  '  c — )/mr  ‘  (6) 

The  canonical  angular  momentum  P0  is  a  constant  of  the  motion  and  the 
magnetic  vector  potential  has  two  components,  i.e., 

A0(r)  =  A*xt(r)  +  A*elf(r)  . 

The  envelop  (boundary)  equation  for  a  charged  particle  in  the  region 

self 

-  6j  <  z  <  0  can  be  easily  determined  when  <f>(r)  =  A0  (r)  =  0.  Substi¬ 
tuting  Eqs.  (4)  and  (6)  into  Eq.  (5)  and  taking  Pp  =  Pz  =  0,  the  boundary 
is  g i ven  by 1 6 
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R±=  2(8l/*)r1f-3±[l  +  (2/6,)^]^},  (7) 

where  *  Pq^/H,  ft  =  qB^mc,  rx  =  (H/2m^ ) ^ ,  and  PQ  >  0  for  q  >  0. 
Equation  (7)  shows  that  crossing  of  the  upper  (+  sign)  and  lower  (-  sign) 
boundaries  in  the  region  -  6X<  z  <  0  occurs  when  1  +  (z/6^  )ai  =  0.  There¬ 
fore,  the  boundary  remains  open  and  thus  a  rotating  particle  can  pass 
through  the  cusp  if 

1  >  0.  (6) 

According  to  Eq.  (8),  a  ring  that  reaches  the  entrance  of  the  cusp  (z  = 

)  with  small  self  fields  and  small  radial  velocity,  it  passes  through 
provided  that 


(m/2)  (v®^  +v^)>  R1Q1(mv^  -  ql^  /2)  , 


where  is  the  major  radius  of  the  ring  at  z 
Eq-  (9)  gives 


V  Since  vei  = 


v2l  >  0, 


which  is  not  an  unexpected  result.  For  v  »  0,  i.e.,  for  a^1  »  1,  Eq. 
(7)  becomes 


R  (z)  = 


'  (1+61r1/|z|  . 


The  discussion  on  the  transmission  of  the  rotating  beam  through  the 
half  cups  is  based  on  the  assumption  that  the  self  magnetic  field  of  the 
beam  is  very  small  in  comparison  to  the  external  field.  However,  the  ex¬ 
istence  of  equilibrium  in  the  region  that  is  located  to  the  right  of  the 
half  cusp  (z  >  0)  requires  the  presence  of  self  fields.  Therefore,  if  the 
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self  magnetic  field  is  zero  at  the  exit  of  the  cusp  means  muct  be  found 
that  inhibit  the  flow  of  electrons  in  both  the  azimuthal  and  axial  di¬ 
rections.  Since  in  multiple  elastic  scattering  the  square  of  the  mean 
scattering  angle17’18  of  charged  particles  injected  with  the  same  velo¬ 
city  into  neutral  gas  is  proportional  to  R2  .cnRT  1 ,  where  R  (=  q2/mc2)  is 

o  u  o 

the  charged  particle  classical  radius,  the  electron  current  can  be  ap¬ 
preciably  inhibited  by  passing  the  beam  through  a  neutral  gas  of  approp¬ 
riate  pressure. 

The  dotted  line  in  Fig.  2  illustrates  qualitatively  how  the  field  in 
a  linear  half  cusp  is  modified  by  the  self  field  of  a  long  rotating  beam. 
Clearly,  the  half  cusp  involves  into  a  full  cusp. 

The  dynamics  of  a  neutralized,  rotating  ion  beam  passing  through  a 
cusp  are  very  complex.  There  is  experimental  evidence19  that  a  fraction 
of  the  light  electrons  is  tied  to  the  field  lines  and  does  not  follow  the 
ions  across  the  cusp.  The  transmission  of  a  rotating  beam  through  a  cusp 
is  presently  under  investigation  with  a  computer  simulation  code. 

For  pellet  fusion,20’21  a  full  cusp  may  be  more  desirable  than  a  half 
cusp.  Again,  it  is  assumed  that  for  z  >  0  the  cusp  is  described  by 

A0  =  rzB3/268,  z  >  0  (12) 

where  6  is  the  transition  width  and  El  is  the  axial  magnetic  field  B  = 

3  3  Z 

\  3 

—  (rAQ)  at  the  exit  of  the  cusp,  i.e.,  of  z  =  6  .  Similarly  to  the 

r  d  r  o  2 

z  <  0  case,  the  boundary  equation  is 

\  =  2r2 (fig/z)  |  T  1  +  [1  +  (z/ss)oa]2  }  ,  z  >  0  (13) 

where  a ^  =  otj  B3/B1  and  r2a2  =  P0/(2mH)^.  Figure  3  shows  the  boundaries 
of  the  full  cusp  for  a1  -  0.5  and  a2  =  0.1,  i.e.,  for  B3  =  B1/5.  Higher 
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values  of  B3  may  be  more  desirable.  A  pellet  that  is  located  near  the 
exit  of  the  cusp  is  irradiated  not  only  from  the  front  but  also  from  the 
back  because  of  the  reflected  ions.  The  range  „f  100  M eV  protons  in  gold 
is  about  0.8  cm  and  therefore  the  total  diameter  of  the  target  (DT  and  Au) 
should  be  ~  1.7  cm.  Thus,  heavier  ions  are  more  suitable.  The  size 
of  the  pellet  can  be  made  compatible  with  the  opening  of  the  boundary 
by  adjusting  the  value  of  B..  . 

IV.  Equi librium 

When  a  rotating  beam  that  is  immersed  in  a  magnetic  field  is  extracted 
to  a  free  magnetic  field  region,  as  in  Fig.  1,  the  canonical  angular  momen¬ 
tum  can  only  be  conserved  if  r  increases,  i.e.,  when  an  expansion  of  the 
beam  occurs.  For  intense  rotating  beams  that  are  not  current  neutralized 
the  vector  potential  Ag(r)  has  two  components, 
i  .e.  , 

Ag (r )  =  AgXt (r)  +  A*eIf(r), 

where  Ag  (r)  is  due  to  the  externally  applied  field  and  Ag  (r)  is  due 
to  the  azimuthal  current  of  the  bear.,  Therefore,  Pg  can  be  conserved  with¬ 
out  an  appreciable  increase  of  r,  even  in  the  absence  of  an  external  field, 

provided  that  A.  (r)  ^  0.  However,  conservation  of  PQ  does  not  warranty 

0  0 

the  existence  of  equilibrium.  For  the  equilibrium  to  exist,  a  negative 
force  is  needed  which  can  be  conveniently  provided  by  the  self  Bg  field. 

For  a  solid  rotating  beam  propagating  in  the  absence  of  an  external 
magnetic  field,  the  J  B„  force  balances  the  centrifugal  and  the  JQB  forces. 

Z  u  u  Z 

This  equi 1 ibri um  was  studied  initially  by  Yoshikawa,22  who  used  a  cold  fluid 
model  (VP  =  0)  and  assumed  b  -*■  «>  (no  surrounding  conducting  wall). 
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However,  only  hollow  beams  can  be  efficiently  transformed  with  mag 


netic  cusps.  For  this  reason  we  study  in  this  section  the  equilibrium 
of  a  hollow  beam.  Any  realistic  model  dealing  with  finite  thickness 
hollow  beams  should  include  the  pressure  gradient  term.  Hot  beams  (p  ?  0) 
can  be  conveniently  studied  using  the  Vlasov-Maxwel 1  approach. 

The  present  calculation  is  carried  within  the  Vlasov-Maxwel 1  frame¬ 
work  and  indicates  the  existence  of  equilibrium  for  a  hollow  rotating 
beam,  when  the  external  magnetic  field  is  zero  and  the  radius  of  the 
conducting  guide  tube  b  -»  °°. 

In  addition  to  these  theoretical  predictions,  recent  experiments23 
at  NRL  have  demonstrated  the  existence  of  equilibrium  for  a  hollow,  over- 
dense,  rotating  electron  beam  that  propagates  inside  a  conducting  guide 
tube  in  the  absence  of  an  external  magnetic  field.  Presently,  experiments 
are  planned  to  test  the  existence  of  equilibrium  when  the  conducting  tube 
that  surrounds  the  beam  is  removed. 

The  equilibrium  configuration  is  shown  in  Fig.  4.  It  consists  of  a 
space  charge  neutralized,  hollow,  infinitely  long  rotating  ion  beam  that 
is  immersed  in  a  uniform  magnetic  field  Bq.  The  beam  is  surrounded  by  a 
conducting  cylindrical  liner  of  radius  b  and  propagates  with  a  velocity 
vz(r)  along  the  axial  direction.  The  following  simplifying  assumptions 
are  made  in  describing  the  equilibrium  of  the  propagating  ion  beam: 

(i)  Equilibrium  properties  are  azimuthal ly  symmetric  and  indepen¬ 
dent  of  z-coordi nate ,  and 

(ii)  The  equilibrium  electric  field  of  the  beam  is  neutralized  by 
an  electron  background  and  the  current  carried  by  the  electrons  is 
equal  to  zero,  i.e.,  the  self  magnetic  field  is  generated  entirely 


by  the  ions. 

Initially,  the  equilibrium  properties  of  the  beam  are  computed  for 

finite  B  and  b.  Subsequently,  the  limited  case  of  B  =0  and  b  00  i  s 
o  o 

considered. 

The  equilibrium  is  described  by  the  distribution  function 


where  to,  1^,  m,  n,  c,  m  and  Pzq  are  constants,  A^(r)  is  the  axial  component 
of  the  vector  potential, 

H  =  (m/2)(vzr  +  v>  +  v*),  (15) 

Pg  =  mrVg  +  (q/c)rAg(r),  (16) 

and 


Pz  =  mvz  +  07) 

The  argument  of  the  first  delta  function  in  Eq.  ( 1 ^ )  may  be  expressed  as 

H  -  ‘"Pg  -  K  =  (m/2)[(v0-r<u)2  +  v2]  +  UQ(r),  (18) 

where 

U  (r)  =  (m/2)[P  /m  -  (q/mc)A  (r)]2 
o  zo  z 

-  (q/c)rA  (r)au  -  mr2  oo2/2  -  k.  <  0.  (19) 

0 


The  density  profile  corresponding  to  the  distribution  of  Eq.  (14)  is 


ajo) 


n(r)  =  f'f° (H-ujP  ,  v  )dav  =  n  ,  a  <  r  <  an  , 

0  z 

where  ^  3nd  a.,  are  the  roots  of  Uo(r)  defined  in  (19). 

The  azimuthal  current  density  is  given  by 

JQ(r)  =  qnV  (r)  ,  a  «  r  5;  a  , ,  (21) 

where  V^(r)  is  the  mean  azimuthal  velocity  of  the  rotating  beam,  i.e., 

V  (r)  =  <v  (r)>  =  (1/n)  [V  f° (H-wP  ,  v  )d  v 
0  0  J  0  0  z 

=  Wr  ,  a,  <  r  <;  a=  .  (22) 

Equation  (22)  shows  that  the  ion  beam  rotates  with  a  constant  angular 
velocity  tv  around  its  axis  of  symmetry,  i.e.,  the  beam  is  in  a  rigid- 
rotor  equilibrium. 

The  azimuthal  component  of  the  magnetic  vector  potential  is  deter¬ 
mined  from 

d7  r  d7  £rVr)]  =  *  (Wc)qnujr,  <23) 

and  the  axial  component  from 

Id  /  dAz(l ^  j  =  -  (4n/c)qn[Pzo/m  -  (q/mc)Az(r)] .  (24) 

r  dr  l  dr  J 

The  solution  of  Eq.  (23)  subject  to  the  boundary  conditions  that 

b 

A  end  B  are  continuous  at  a.  and  a_  and  [2^15  rdr  =  (flux  con- 

9  z  ~  a  z  o 

servation)  is 

IBz(a1  )r/2  ,  0  ^  r  ^  a, 

Bz(a1  )r/2  -  Trnqu)(r2-a^)2/2cr  ,  a1<  r  a2  (25) 

Bz(a1)r/2  -  nnq<u[  (a2-a^) 

X  (2r2-a^ -a2)]/2cr ,  ap  <;  r  <  b, 
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where  B  (ax)  is  the  axial  magnetic  field  in  the  inner  surface  of  the  beam. 
The  axial  magnetic  field  Bz(r)  =  "jr  gp(Agr)  *s  computed  from  Eq.  (25) 


where  I  =  3^(a2  -  a2)  is  the  azimuthal  current  per  unit  length.  For 
q  >  0,  1^  has  the  same  sign  with  oo,  which  is,  either  negative  or  positive. 

Although  the  expressions  describing  the  density,  azimuthal  velocity, 
Ag(r)  and  Bz(r)  appear  to  be  identical  to  those  of  reference  2k,  they  are 
defined  in  a  different  region  because  for  the  distribution  function  of 
Eq.  (14),  ax  and  a^  are  functions  of  A^lr). 

The  axial  component  of  the  vector  potential  A^(r)  is  determined  from 
Eq.  (2*0  and  is  given  by 


Az(r)  = 


_  ■  _  “i 


(cP^/qHl-U/A*)  [l0(r/A*)K1(a1/A*)  +  K0(r/A*)I1(a1/A*)]}f  ai<r<a2, 
(2I/c)£n(aa/r)  +( cPZQ/q ) {1  - ( a  ,/a*)  [lQ ( a, / A*) Kq( /a*)  +  <27) 

K^a./A^IJa^A*)]} 


where  I  and  K  are  modified  Bessel  functions  of  order  n,  A*  =  c/u) 
n  n  p 

and  w2  =  ATrnq2/m.  Equation  (27)  has  been  derived  assuming  that  A  and  B, 


are  continuous  at  a,  and  a_  and  A  (a  )  *  BQ (a,  )  =  0.  The  azimuthal  magnetic 

A  3  2  1  U 

field  Bq  =  -dA  /dr  computed  from  Eq.  (27)  is 
o  z 
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(28) 


0,  Osrs  a, 

BQ(r)  =  <^qn.lP,0/«»c)[Kl(al/x*)Ix(r/x*)-J1(a1/x*)KJ(r/x*)],asr«3 

21/ cr ,  a,  5  r  $  b, 

where  the  current  of  the  beam  I  can  be  expressed  as 

I  =  (2nqnP2oa1aa/m)[K1(a1/\*)I1  (a3/x*)  -  \  (a1  /x*)^  (r/x*)]  .  (29) 

For  a^/X*  »  1,  Eqs.  (27),  (28)  and  (29)  become 

0,  0  s  r  ^  a1 

Az(r)  =  (cPzo/q)fl  -  (a1/r)^  cosh[(r-ai  )/x*]} ,  aa  s  r  <;  aa  (30) 

(2I/c)tn(a2/r)  +  (cPzQ/q)[l  -  (a1/a2)|  cosh[ (a.,^  )/\*]\  &  ^r<b 

0,  0  ^  r  s  a, 

BQ(r)  =  (4nqnx*Pzo/mc) (a1 /r)^  sinh[ (r-a1 )/x*] ,  a1  <  r  s  aa  (31) 

2I/cr,  aa  s  r  s  b, 

and 

I  =  (2nqnPzoX*/m)  (a1aa)^sinh  [  (a.^  )/x*]  •  (32) 

|  Defining  the  Alfven  limiting25  current  as  the  beam  current  at  which 
the  Larmor  radius  of  a  charged  particle  in  the  maximum  self  magnetic  field 
is  equal  to  one_ha!f  of  the  beam  radius,  we  obtain 

W _ s 

(q/mc) (2IA/aac)  c  *  (33) 

For  ax/X*  »  1,  Eqs.  (30)  and  (33)  give 

lA  =  c2(P2o/q)(ai/a2^  cosh[(aa  -  a1)/x*].  (34) 

Combining  Eqs.  (32)  and  (34),  we  get 

I  as  * 

~  =  2X*  tanh  £(aa  -  a1>/X  ]» 

A 
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which  is  reduced  to 


J_  =  a.^/2\* ,  for  (aa  -  ax  )/A*  »  1  (3k) 

JA 


The  conclusion  from  Eq.  (3**)  is  that  the  beam  current  I  can  exceed 
considerably  the  Alfven  current  1^  when  the  collisionless  skin  depth  X* 
is  much  smaller  than  the  radius26  and  the  thickness  of  the  beam. 

The  function  -UQ(r)  >  0  has  a  maximum,  which  can  be  determined  from 


Using  Eq.  (19),  Eq.  (35)  can  be  written  as 


"  (?)Vz(po)B0(Po)  +  (c)u,poBz(po)  +nuj2P0  =  0 


or 


where  ft  (p  )  =  qB  (p  )/mc  and  ftn  =  qBn(p  )/mc.  Since  the  radical  in  Eq. 
zo  zo  o  tto 

(36)  is  greater  than  unity,  the  ratio  to/ftz(pQ)  can  be  either  positive 
or  negative.  Which  mode  of  rotation  is  present  depends  upon  the  initial 
preparation  of  the  system. 

When  vz  =  const.  Eq.  (36)  gives 

o>  =*  -  ^z(po)  . 

as  expected.27*20  It  can  be  shown29  that  the  pressure  P  in  the  beam  is 
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(37) 


P  =  -nU  (r) . 
o 


Using  the  1-fluid  equation  with  ^  =  0,  we  obtain 


nm(V'V)V  =  nq  -■  *  ^  -  Vp, 


or 


V2 

-  4  ™  ‘  (j6b2  '  Jz6e>/C  •  37-  <38) 


The  balance  of  forces  becomes  apparent  when  Eqs.  (19)  and  (37)  are 
inserted  in  Eq.  (38). 

The  inner  a  and  outer  a  radii  of  the  beam  are  the  roots  of  U  (r) 
13  o 

function  and  therefore  U  (a,  )  =  U  (a  )  =  0.  These  radii  are  in  general 

o  A  o  a 

complicated  functions  of  several  parameters.  Of  particular  interest  in  the 

present  work  is  the  limited  case  B  =  0  and  b  00 

o 

Defining  the  layer  strength  parameter  as 

V  =  Ttnq2  (a2  -  a2)/mc2,  (39) 


the  constant  B  (a1 )  in  Eq.  (25)  can  be  written  as 


Bz  (ax  )  =  2u)mcv/q . 


(**0) 


and  the  azimuthal  component  of  the  vector  potential  Ag(r)  in  the  region 
\  —  r  —  aa  as 

]•  (4l) 


A0(r)  =  (muurc/ q)i  1  - 
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Substituting  Eqs.  (41)  and  (17)  in  Eq.  (19),  it  is  obtained 


UQ(r)  =  U(r)  +  mv2(r)/2, 


where 


U(r)  = 


\*u?ut 


+  i  >2*,r2  + 


1 


(a;-af)J  -2(a^-af) 


<  -  \  I,  ^ 


and 


vz(r)  =  v2(ai  )(ai/X*)[1c(r/X*)K1  (a1/x*)  +  ^(r/x*)!,  (a,  /X*)]  (44) 


The  constant  vz  (a^  )  =  Pzq  /m  is  the  velocity  of  charged  particles  in 
the  inner  radius  of  the  beam. 

The  functions  U(r)  and  -mv2(r)/2  are  sketched  in  Fig.  6.  The  two 
roots  a^;  and  a^’  of  U(r)  are  given  by 

-  o  ♦  *,)<<£  -  *  *1 

*  Ctt  +  2^’  <a='“i,)  +  *>]  '  [“>' 

The  maximum  value  of  Ll(r)  occurs  at 

p2  =  a^  +  (a^  -  a^)/2v  =  a^  +  (2TrnRQ)  1 


and  i  s 


U(r> ) 


t!KU2(2V  +  1) 


=  [ af(2v  - 


1)  +  a2(2v+l)J/8v+k1/iiiD2(2v+l),  (45) 
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- - »-■  '  -  r  i  IT— i 


where  Rq  *  q2/mc2  is  the  charged  particle  classical  radius.  Necessary 
conditions  for  the  equilibrium  to  exist  are 

U(p)  <  0 

and 

|  U ( p>  |  >  |  v|(aa)  .  (46) 

Substituting  Eqs.  (44)  and  (45)  in  Eq.  (46),  it  is  obtained 

4vv*(«  )((a  /\*Tl  (a  /X*)K  (a  /X*)+K  (a  /\*)I  /x*f!*  -  1}  (47) 

(2^1)3  >  _ _ _ 1 _ -  °  5  11  °  S  1  -* 

<*>2(af-af) 

The  radii  of  the  beam  a i  and  ag  are  the  roots  of  UQ(r)  and  are  deter¬ 
mined  from  Eq.  (42).  These  two  radii  are  given  by 

ai  =  ^Pzo  "  2mki)/n,2(*)2(1  +  2v^’  ^ 

and 

a,2  (v  +  1 )  +  va2  +  2k^  /mu>2  * 

Tv^aj  )/a)2J(a1/X*)^”lo(a,/x*)K1  (a1 /x*)  +  K,^/**)1!  (a1/x*)T-  (*9) 

When  ax/X*  »  1  and  (a3  -  aJ/X*  »  1,  Eq.  (49)  becomes 

a  CD2  _ 

(a2-a^»  ~  (X*/2)^n  {— - I  (^Da2  +  vaf  +  2\  !m?  . 

a1v2  (a,)  - 

In  order  to  have  a3  ■  a^  ^  0, 

jjo2/v2(a1  )"[<v  +  1  )«f  +  vaf  +  2kj  /nW2J  ;>  ^  /a2  . 
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1 


(b)  Infinitesimally  thin  beam 

The  various  expressions  given  above  are  simplified  considerably  when 
a  ~  a  ,  i.e.,  for  very  thin  beams. 

In  the  thin  beam  approximation  the  radial  velocity  of  particles  is 
approximately  zero  and  then 

fv2(r)  H  -  jrav2  =  H  -  (m/2)£p0/mr  -  (q/mc)Aj  (50) 

Using  Eq.  (50)  and  the  relation  H  -  u>Pg  =  1^,  Eq.  (42)  gives 

a*  a  a2  =  Pg/mw(1  +  v).  (51) 

Equation  (51)  can  be  easily  derived  from  the  conservation  of  PQ  = 

0 

mr  e  +  qAQr/c.  Substituting  6  -  co  and  using  AQ  =  mvcorc/q  from  Eq.  (4l), 
we  get  Eq.  (51).  For  a^/X*  »  1,  Eq.  (47)  becomes 

(2v  +  l)2  >  4nv^(al  )nRQ{(al/a3)  cosh2[  (a2-a1  )/\*]  -  1)  of2  . 

Since 

I  =  2nqnvz<a1  ^aj  asX  sinh[a5  -  a1)/\*],  in  the  thin  beam  approximation 
Eq.  (47)  becomes 

(2v  +  Djr"  -  (1/a))'  (52) 

o 

The  above  equation  shows  that  in  thin  beams  the  current  T  can  be 
large  provided  that  the  frequency  of  rotation  is  high. 

Equation  (52)  can  be  written  in  a  different  form.  Since  I  -  qvV  /R 
and  wa^  =  vQ,  Eq.  (52)  gives 
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Since  for  infinitesimally  thin  beam  ft2(po)  =  wv  and  ftg(pQ)  =  qx/mc2^, 
the  equation  that  precedes  Eq.  (36)  gives 


V6  _  q  x _ 

Vz  mc2ai  (v  +  1  )oj 


(54) 


or  after  substituting  for  X 


ve  _  /  v  \j  (55) 

vz  '  V  1  +  V  ) 


The  inequality  of  Eq.  (53)  is  always  satisfied  when  the  equilibrium 
condition  of  Eq.  (55)  is  satisfied. 

When  the  distance  between  the  accelerator  and  the  target  is  much 
shorter  than  the  beam  length,  a  steady  state  could  be  established.  Under 
these  conditions,  using  Eq.  (15)  with  vr  =  o  and  Eqs.  (51),  (54)  and 
Vg  =  we  obtain 


(v  +  l)2v 
(2v  +  1) 


2maaH 


i  P 


(56) 


where  Hp  is  the  hamiltonian  at  the  peak  of  the  compression.  Before  the 
establishment  of  the  steady  state,  a  fraction  of  the  beam  energy  is  in¬ 
vested  to  build  up  the  magnetic  energy.  If  L  is  the  beam  inductance  per 
unit  length,  then  Eq.  (56)  becomes 


p2 

(V+l)g  V  _  =  fl 

[(2V+1)  +  Lvcg  (v+-l)]  2maf  Hp 


(57) 
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The  beam  inductance  L  has  two  terms,  one  associated  with  the  rotation 
and  the  other  with  the  translation  of  the  beam.  In  general,  L  is  a  func¬ 
tion  of  v.  For  an  infinitesimally  thin  beam,  of  length  1,  the  inductance 
per  unit  length  is  given  by 

L  -  c-2[21n(2l/ai)-l  +  ] 

provided  1  »  a  . 

i 

(c)  An  Example 

Combining  Eqs.  (51),  (55)  and  (56),  it  is  obtained 


V2  (\>+l ) 

(2\>+T)  ’ 


(58) 


where  Vqq  is  the  azimuthal  velocity  of  ions  at  the  peak  of  the  compres¬ 
sion  (H  =  imvQ2).  Equation  (58)  determines  \>  when  the  value  of  I  and 

p  DO 

Vqq  are  specified.  For  I  =  107A  and  =  200  MeV,  Eq.  (58)  gives  V  =  0.8. 

The  value  of  V  is  determined  from  V  =  IR  /qv  and  is  1.19X1010  cm/sec. 
z  z  o 

Substituting  this  value  in  Eq.  (55),  we  get  Vg  =  0.8  X  1010  cm/sec. 

The  beam  radius  is  determined  from  Eq.  (51).  Using  the  value  of 
PQ  that  corresponds  to  the  values  of  parameters  of  Table  1,  Eq.  (51)  gives 
a^  0.5^  cm.  Since  for  this  particular  case  the  length  of  the  extracted 
beam  is  only  38  cm,  the  ommission  of  inductive  effects  is  justified 
only  if  the  target  is  located  a  few  cm  from  the  exit  of  the  cusp. 

In  the  above  example  the  beam  has  been  assumed  to  be  infinitesimally 
thin.  If  the  thickness  of  the  beam  is  finite,  the  procedure  of  deter¬ 
mining  the  various  parameters  of  the  beam  is  more  complicated.  Finally, 
it  should  be  emphasized  that  predictions  that  are  based  on  infinitesimally 
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thin  beams  must  be  treated  with  care,  because  such  beams  are  likely  to  be 
unstable. 

V.  Summary 

A  scheme  is  proposed  for  the  generation  of  high  energy,  high  current 
ion  pulses  by  magnetic  compression  of  ion  rings.  An  important  feature 
of  the  present  scheme  is  the  extraction  and  unwinding  of  the  ring  after 
compression.  The  propagation  of  the  extracted  pulse  from  the  compression 
region  to  the  target  does  not  require  the  presence  of  metallic  boundaries 
or  the  application  of  an  external  magnetic  field.  Therefore,  such  a 
scheme  is  in  particular  useful  for  targets  that  are  separated  from  the 
compressor  by  a  large  distance.  It  appears  that  the  present  technique 
can  lead  to  the  generation  of  hundred  of  MeV  ion  pulses  without  the  need 
for  further  development  in  the  existing  pulse  power  technology.  The 
capabilities  of  the  proposed  scheme  appear  to  be  far  beyond  the  present 
requirements  for  pellet  fusion.  However,  the  successful  evolution  of  the 
proposed  approach  to  a  practical  device  rests  very  heavily  or,  several 
factors,  including  the  formation  and  stability  of  ion  rings,  the  stability 
and  finite  length  of  the  hollow  beam  and  the  ability  to  inhibit,  at  least 
partially,  the  electron  return  current  in  the  propagating  ion  beam. 
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Table  1 


Source 

Average  hollow  beam  radius  at  the  anode  (cm) 
Magnetic  Field  (kG) 

Pulse  Duration  (nsec) 


26 

2.3 

70 


Compression  Region 


Proton  energy  (MeV) 

Magnetic  field  at  M.P.  (kG) 
Proton  current  (MA) 

Number  of  protons 
Ring  kinetic  energy  (MJ) 


Before 

Compression 

2 

18.2 

0.92 

2  X  1017 
0.064 


After 

Compress i on 

100 

910 

46 

2  X  1017 
3.2 
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Fig.  3  —  Boundary  of  the  allowed  region  (hatched)  of  a  full  cusp  when 
oj  =  0.5  and  a 2  =  0.1 


nm(V*V)V  =  -  VP 

nmv|/r  +  ^[jeBz-JzBe]--^=0 

Fig.  5  —  Illustration  of  various  forces  acting  on  the  beam 


